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( We consider spontaneous magnetic flux in a doubly connected s- wave/normal metal/d-wave 

superconductor (SND) Josephson junction. The flux magnitude is not quantized and is determined 
by the parameters of the system. In an isolated system, quantum fluctuations of the superconducting 
phase take place, leading to quantum flux noise. The possibilities of experimental observations of 

■ the effect are discussed. 

Mesoscopic SNS systems containing superconductors with different pairing symmetry promise a whole new range 
of phenomena, such as half-periodic Josephson effect [Q and spontaneous currents parallel to the interface . In the 
latter paper, it was shown that due to local violation of T-symmetry, spontaneous currents should flow parallel to the 

■ SN interface in the normal part of the system; the authors noted that such a current is difficult to observe directly, 
since it will average to zero over short distances. 

This difficulty does not appear in a doubly connected SND system. There the spontaneous current would give 



rise to an equilibrium magnetic flux. Its magnitude is less than <J> and is determined only by the parameters of the 
system. It should not be mixed up with the fractional flux that is predicted in a contour with a T-symmetry- violating 
junction (^j-Q]. The latter is created by the current flow across the junction; the former is directly created by the 



spontaneous current flowing parallel to the SN interface. It is therefore an analog of persistent currents in normal 
mesoscopic rings ^] , when the T-symmetry is broken by an external magnetic flux. 

We consider here an annular SND junction (Fig.|l|). This idealised structure consists of the inner (s-wave) supercon- 
' ducting contact S, conducting ring N, and the outer (d-wave) superconducting ring D. The splits in S and D ensure 
that the screening supercurrents in these electrodes do not mask the effects of spontaneous current in the normal part 
of the system, N. For the sake of simplicity, we assume that the orientation of the d-wave order parameter positive 
and negative lobes with respect to the radius vector is the same along the d-wave ring. 

We calculate Josephson current in the system, following Q . It is carried through the normal region by two sets of 
Andreev levels, coupled to positive and negative lobes of the d-wave order parameter (chosen to be real; the phase 
difference across the contact is thus (f>, where Aie 1 ^ is the order parameter of the s-wave superconductor). We are 
considering a "small" junction limit, where the influence on the Josephson current of its own magnetic field can be 
' neglected. This imposes a condition (see Fig.l) 



o 



2irR < Xj, (I) 



where Aj = \/frc 2 /87rej c A is the Josephson penetration length ||; j c is the critical Josephson current density, 
and A = W + + Xl,2 is the integral London penetration length, Al i0 being the penetration length in the ath 
superconductor. In case of SNS junction 



o 

^ I where Xf is Fermi wavelength in the normal part of the system, and a/ s = e 2 /Tic ~ 1/137. 

Using the set of cylindric coordinates (r, a, z), we write the Bogoliubov-de Gennes equations in the normal region 
i for a mode with vertical momentum k z as 

f-^F^ 2 -^ ° i 2 \ f u(r,a;k z ) \_ £ fu(r,a;k z )\ ^ 



2^tV 2 + fi±J \ v{r,a;k z ) J \y(r,a;k z ) 

where /Ltj_ = (k F — k z )/2m* = fc^/2m*, and m* is the effective mass of the electron. Expanding u,v over asymuthal 
modes, u(r, a) = J2 n Un( r ) ema /Vr, v(r, a) = J2 n Un( r ) ema /Vr, we find: 

rf n 2 — 1 I A 

- ^2 U t(r)+—^Ut(r) = (±2m*E+ki)Ut(r). (4) 
The Andreev levels are determined from the quantization condition |jj 
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pdr= drp(r,E)+ / dr p(r,-E) = 2nm + n ± (<j) + S n , k J, (5) 
Jr Jr+w 

where m is integer, and j?(r, E) = (kj_ + (n 2 — l/4)/r 2 + 2m* E) 1 / 2 . The intrinsic phase shift 5 n>k in mode (n,k z ) is 
zero if the corresponding lobe of d-wave order parameter carries positive sign, and n otherwise!. The radial and 
asymuthal currents in a n, fc^-mode are now easily found (cf. ^,||) as: 



Here ^ is the partial radial (asymuthal) current carried by Andreev levels with energies BsiO: 



■r , t, s - a ee(n, k z )n f 

Jn,k, ~ ee ( n ^z); 3n,k, ~ 27r(n 2 - 1/4)1/2 I arCtan 



(k ± {R + W) ) 
n 2 - 1/4 
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— 1 



— arctan 



(k x R) 
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n 2 - 1/4 



- 1 



(7) 



and 



k 2 _/2m* 

£(n ' *' = [k\(R + W) 2 - (n 2 - 1/4)] 1/2 - [jfc2 i?2 _ („2 _ 1/4)] 1/2 ( ) 

plays the role of the interlevel spacing v F /2W of a planar contact. The currents depend on the phase through 
F( * 1 , * (-l) m+1 sinm(0 + S n . kz ) ne- m ' r « n ^ 

r[ (3e(n,k z yre(n,k z y 9 + ° n - kJ n sinh /3c(»,fc a ) ' W 

where r 3> e(n, fc 2 ) _1 is the elastic scattering time. In the limit of zero temperature and no scattering F(0,0;4>) 
becomes the 27r-periodic sawtooth of unit amplitude. 
First we find the Josephson (radial) current: 

+ 

l](4>) = J2 Jr ( n > + E Jr ( n > ^ + ""J- (10) 

n,k z n,k z 

The "±"-sums are taken over zero- (5 n k z = 0), and 7r-levels (5 n ,k z — t), coupled to positive and negative lobes of the 
d-wave order parameter respectively Q . At zero temperature and without elastic scattering this reduces to 

Ij{4>) = ^-I F{0, 0; <j>) + l^-I F(0, 0;<t>+ n), (11) 

where Iq = N±ee is the critical current in the system, N± ~ A/ \p is the number of quantum conducting channels 
through the normal part of the system (of crosssection area A), and e = (e(n, k z )) ~ Vp /2W is the average interlevel 
spacing. The imbalance coefficient Z (\Z\ < 1) depends on the orientation of d-wave superconductor: 

z = lJ2 ^ - E < n > k ^ | / 1 E e ( n ' fc *) + E e ( n > | . ( 12 ) 

and determines the equilibrium phase difference across the junction, 0o = ±-^^-tt. The Josephson energy of the 
system is thus 

17(0 = 1 / = (-7T < < 7T). (13) 

In equilibrium, the total Josephson current is zero, because contributions from zero- and 7r-levels cancel. On the 
contrary, their contributions to the asymuthal current add up and lead to finite spontaneous current. In the situation 
of Figjl], due to partial cancellation of terms with positive and negative n, it equals 

Ek.T^n-Ce,k M )+iW\n\,k,-,<f>)- J a (\n\,k z ;cf> + n)) if n+(0, k z ) > n -(0, k z ); 
l s (<f>) = {0 if n+(0,fc») =n-(6,k z ); (14) 

- Efe, E^l%] kz)+1 (Ja(\n\, h; 4>) - J a (\n\,k z ;cf> + it)) if n+(9, k z ) < n~(0, k z ). 
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Here ([x] means integral part of a;) 



n+(6,k z 



tan6» v /(fc ± i? 2 ) 2 + 1/4 
^l + tan 2 e»(fc ± i? 2 ) 2 



n-(0,k z ) 



cot6'V(fc_ L fl 2 ) 2 + l/4 
v/1 + cot 2 6(k ± R 2 ) 2 



(15) 



In the limit /?, r, i?, R + IF — * oo we return to the result fa] for a planar junction 



/ s (0) = ^^^-sin(0-^)(F(O,O;0)-F(O,O;0 + ^)); < < tt/2. (16) 

In the annular junction spontaneous current produces a spontaneous magnetic flux, $ s = —I s , where C ~ 7ri? is the 
self-inductance of the system. Flux magnitude can be estimated as 

$ (17) 

The phase difference <j> = \i — X2 is canonically conjugate to the operator of difference of Cooper pairs' numbers in 
the superconducting banks, An = fi\ — fi2 J§]. Therefore in the presence of charging energy due to electron transfer 
between the banks in an isolated junction (with classical capacitance C), <fi ceases to be a constant of motion. Its 
evolution can be mapped on the motion of a quantum particle in a one-dimensional ring — ir < (f> < tt (see e.g. || and 
references therein), with the Hamiltonian 

d . ft 2 A d N2 TT , „ 1# Cfi 2 fi 2 
o<p 2Mq i 0(p 16e z »£q 

The first term describes the charging energy. The effective potential U{4>) (Fig.||) corresponds to the Josephson 
energy: 

m = - M2 = ^ m - 00)2 = ^ m - M2 > (i9) 



where e was introduced earlier, and ujq = y 1 [j^j ' = y ~^sN± ' ipr is the frequency of small oscillations in a 

separate well. The dynamics of the system is thus determined by four parameters: sq, e, N±, and 4>q. We will consider 
the case \4>q\ <C tt (or \tt — 4>o\ <C 7r, in which case we should take < cj) < 2ir). Then in the limit Tiluq, ksT <C U(Q), 
we can use the two- well approximation fici[| , and the transition rate between the wells will be given by T — va + vt, 

where the rate for thermally activated transitions va ~ uj e and for quantum tunneling ~ <^o e ni "° • 

The phase fluctuations create finite voltage, V = w-4>„ and therefore normal current GV flows in the normal region 
(its conductance G ~ N±-^-). The corresponding dissipative function and the decay decrement are 

^4^4^4 G (A)> ; 1= J^& = ®* = ± N J±. (20) 
2 eft 2 2 V 2e / M Q <j>d^> 4e 2 M Q tt ^ ft y ' 

The character of quantum fluctuations of the magnetic flux drastically depends on dissipation. In the limit Hujq <C U(0) 
the conclusions basically reduce to the following The tunneling rate exponent multiplies by ~ (1 + 7/ojq)- 

Therefore if 7 > u>o, tunneling is suppressed. (Strictly speaking, in this case the description in terms of initial 
two- well potential becomes meaningless.) If, on the other hand, 



— = JA^«1, (21) 

UJQ V 6 

quantum tunneling is possible, but quantum coherence between the states in left and right wells is destroyed [Tlfl . 
The problem is simplified due to phase dependence of the spontaneous current (|l6|), which can be approximated by 

$(</>) = $ s sgn0(mod27r). (22) 

This allows to describe the dynamics of <E> in terms of a two-level system p2[ |: state "J," is when (j) is in the left well 
(<f> = — $ s ), state "t" when <fi is in the right well ($ = Q s ). Therefore we can write a set of coupled master equations 
for diagonal terms of density matrix ( O], Ch.7): 



3 



|fPT = r(^-PT) (23) 

Here gives the probability of occupation of right (left) well (spontaneous flux ±$ s ). Eq.(|23|) describes random 
telegraph noise of spontaneous magnetic flux, with autocorrelation function and spectral density p"5| 

43> 2 r 

(*(t)$(Q)) = ^e- 2r l*l; ($ 2 ) w = ~rf^- (24) 

The noise is temperature dependent at T > T* = htoo / fcs , where thermally activated transitions are more prominent 
than tunneling. Below T*, the activated proccssc freeze out, and quantum tunneling remains the sole source of flux 
noise. Magnetic noise intensity detected at a certain frequency J7, as T is changed with temperature, will peak at 

r(T) = n/2. 

Let us make some numerical estimates. Consider the case when the role of the normal part of the system is played 
by 2D electron gas of a GaAs-AlGaAs heterostructure, with \p — 450A and vf — 2 x 10 7 cm/sec. Then the number of 
transverse modes is N± — AirR/Xp, and the condition of "small" contact (||) becomes N± < 2(^y / W/A ss 250-y/W/A. 
It will be satisfied if e.g. N± = 200, (R w 7000 A ), W = 1000 A. The interlevel spacing is then e « 10- 15 erg. The 
corresponding spontaneous flux amplitude $ s 0.02<f>o is within the experimental capabilities. 

The applicability condition of the two-well model imposes the condition eq/£ <§; (N j_</>o) / '1287T 2 , which will be 
satisfied if fa > (1287r 2 /Ar_ L ) 1 / 4 AT7 1/4 » 0.42. 

The weak dissipation condition (|2l|) translates into €q ^> 5x 10~ 18 erg (that is, total capacitance of the system should 

exceed ~3x 10~ 14 F). This yields for hu> <C \J ^-e 2 w e, that is, w <C 10 12 sec -1 . The flux noise characteristic fre- 
quency r is exponentially dependent on (f>^ (through U(0)), and scales with temperature from ut ~ ujq exp[— U(0)/huio] 
at T T* to v T ~ Wo at T»r, and will probably reach a fraction of GHz. 

In conclusion, we have found the Josephson current and spontaneous magnetic flux in an annular SND junction. 
The magnitude of the flux is not quantized; it depends on the parameters of the system and is of the order of few 
percent of <&o for a typical configuration. 

We demonstrated that in an isolated system, due to charging effects, there will be quantum noise of the spontaneous 
magnetic flux, with characteristic frequency strongly dependent on the configuration of the system and its surround- 
ings. It can reach a fraction of GHz. The noise intensity at a fixed frequency will be a non- monotone function of 
temperature, allowing to measure the transition rate T(T) directly and thus separate the contributions from thermally 
activated processes and quantum tunneling. 

We hope that such systems provide experimental observation of a novel dynamical effect in rf-wave superconductivity. 
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FIG. 1. Annular SND junction. The c-axis of the d-wave superconductor is chosen to be parallel to the SD boundary; 9 is 
the angle between the SD boundary and the nodal plane of the d-wave order parameter (0 < 8 < -|). 
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FIG. 2. Effective potential and transition rates for the phase variable in an annular SND junction. 
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